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Abstract. Diluted ferromagnetic III-V semiconductors typically show a high degree of 
compensation. Compensation is connected to the presence of comparable densities of charged 
defects of either sign. This naturally leads to the development of strong correlations between 
defect positions during growth and annealing. We show that these correlations are required 
to understand the experimentally observed transport and magnetic properties as well as the 
persistence of the energy gap upon doping with magnetic ions. 



1. Introduction 

Diluted ferromagnetic III-V semiconductors are promising materials for applications as 
well as interesting from the physics point of view. They might allow the incorporation 
of ferromagnetic elements into semiconductor devices, and thus the integration of data 
processing and magnetic storage on a single chip. It is important to understand the interplay 
between transport properties, magnetic ordering, and the defect configuration, which is 
determined by the growth conditions. We show that correlated positions of defects are 
required for a description consistent with experiments [Jl]]. 

To be specific, we consider GaAs doped with manganese, but the results should apply for 
most III-V compounds. The physics relies on the dual role played by the Mn impurities: They 
carry a local spin due to the half-filled J-shell and dope the system with holes, which mediate 
a ferromagnetic indirect exchange interaction between the spins. A crucial feature of these 
materials is their high degree of compensation, which is presumably due to antisite defects 
(As substituted for Ga). The accordingly small density of holes leads to weak electronic 
screening of the Coulomb interaction between charged defects and between defects and holes. 
This role of compensation has not been taken into account previously. Due to the strong 
Coulomb interactions a random distribution of defects on the cation sublattice is very costly 
in energy. Therefore, defect diffusion 0] leads to a rearrangement of defects in such a way 
that the Coulomb energy is reduced. We study this rearrangement with the help of Monte 
Carlo (MC) simulations and find that strong correlations of defect positions develop. The 
resulting disorder potential acting on the valence-band holes is strongly reduced and most of 
its correlations fall off on the scale of the minimum defect separation. We show that this effect 
has tremendous consequences for transport properties, the integrity of the energy gap, and the 
temperature-dependent magnetisation. In fact, a description in agreement with experiments 
requires the defects to be strongly correlated. 



2. Correlated defects 



In this section, we study how diffusion of defects, i.e., Mn impurities and antisites, changes 
the defect configuration and the disorder potential. This is motivated by experiments [[Zp that 
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show that defect diffusion is rather pronounced at typical growth and annealing temperatures 
of the order of 250°C. We start from a random distribution of Mn impuritites and antisites 
(As) on the cation sublattice, where the density of antisites is determined by charge neutrality 
from the observed |J3p hole concentrations at given Mn doping levels. The Hamiltonian reads 

Z ij tr ij 

where q i are the defect charges and r- • is their separation. Relative to the cation sublattice, 
Mn impurities carry charge q = — 1 whereas antisites carry q = +2. The screening length 
r scr is obtained from nonlinear screening theory [Q. It is much larger than the nearest- 
neighbor separation on the cation sublattice so that it hardly affects the small-scale defect 
correlations relevant here. We perform MC simulations for the Hamiltonian H [J]] employing 
the Metropolis algorithm at the temperature 250°C for systems of 20 x 20 x 20 conventional 
cubic unit cells with periodic boundary conditions, unless stated otherwise. To roughly 
model the dynamics, we only use local MC moves, i.e., exchanges of nearest-neighbor defect 
and/or host atoms. A more realistic modelling of the dynamics involving interstitials and 
vacancies lies beyond the scope of this work. The resulting configurations are quenched at 
low temperatures [Q]. 




Figure 1. Left: Random configuration of defects onalOxlOx 10 lattice for Mn concentration 
x = 0.05 and p = 0.3 holes per Mn. The host lattice is not shown. Small blue circles denote Mn 
impurities, whereas large red ones are As antisites. When two defects sit on nearest-neighbor 
sites, they are connected by a solid line. Right: Equilibrated configuration of defects. It is 
obvious that the defect positions are now strongly correlated. 

Figure [I] (left) shows a random configuration of defects for a Mn concentration of 
x = 0.05 and p = 0.3 holes per Mn, which determines the concentration of antisites (0.0175). 
This should be compared to the right part of Fig. HJ, which shows an equilibrated configuration 
of the same number of defects. It is obvious that the defects have formed small clusters, 
where typically an antisites (charge +2) is surrounded by several Mn impurities (charge —1). 
It is clear that this clustering leads to a strongly reduced impurity potential, since the antisite 
charges are very effectively screened. 

More quantitative information can be obtained from the potential correlation function 

D(r) = (V(r)V(r>)),,-(V) 2 . (2) 
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Figure 2. Potential correlation function D(r) for Mn concentrations = 0.05 and p = 0.3 holes 
per Mn, plotted for various numbers of MC steps. The heavy solid bar denotes the nearest- 
neighbor separation on the cation sublattice. The inset shows the initial rapid decrease of D(r) 
for the equilibrated case, see text. 

Obviously, AV = a/£>(0) is the width of the distribution of V(r). D(r) is plotted in Fig. || 
for x = 0.05 and p = 0.3 for various numbers of MC steps, which serve as a rough meassure 
of annealing time. Defect diffusion strongly reduces D(r) and thus AV. The inset in Fig. || 
shows a blowup of D(r) for the approximately equilibrated state. The initial reduction of D(r) 
to about 10% of D(0) on the scale of the nearest-neighbor separation of about 4.00A(heavy 
bar) is due to the screening of the compensated 70% of Mn impurities by antisites, leaving 
only 30% active. Since D(r) contains the potential squared, this leads to a reduction to about 
9% (yellow region in the inset). The remaining contribution from uncompensated Mn (light 
blue region) decays on the typical length scale of 14.4Aassociated with their density, which 
is 0.3 x 0.05 per cation site. Clearly, ionic screening, i.e., screening by the defects, is nearly 
perfect at 250°C. While the reduction of AV is thus substantial, we find that AV is still not 
small compared to the Fermi energy [[[]] . 

3. Valence-band holes 

In this section we consider the properties of valence-band holes in the impurity potential 
V(r) due to the charged defects. In particular, we are interested in their spectrum and 
localisation properties. We employ the envelope function and parabolic-band approximations 
for the holes and start from the Hamiltonian H = — £-(/? 2 /2m*) V? + V(r i ). For material- 
specific calculations the detailed band structure should be taken into account, e.g., using the 
6-band Kohn-Luttinger Hamiltonian [0]. This should not change the qualitative results but 
is known to increase the mean-field T c [[7]]. The hole Hamiltonian is written in a plane- wave 
basis and diagonalized numerically [JTp . The calculations are done for spin-less holes, since 
the additional disorder introduced by the exchange interaction is found to be much smaller 
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than AV. We obtain the energy spectrum and normalized eigenfunctions \^ n (r). From the 
latter we calculate the inverse participation ratios 

IPRH = l/^|v/, 1 (r)| 4 (3) 

r 

of the states. The main physical content of the IPR is that it scales with system size for 
extended states but essentially remains constant for localized states. The IPR thus allows to 
estimate the position of the mobility edge in the valence band. 
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Figure 3. IPR as a function of electron energy for x = 0.05 and p = 0.3 after various numbers 
of MC steps increasing from the red to the violet curve. The Fermi energies for the random 
and the equilibrated case are also shown. 

Figure [J shows the IPR as a function of electron energy for the same parameters as above, 
x = 0.05 and p = 0.3, after various numbers of MC steps. The plot shows that uncorrected 
defects would lead to the filling-in of the band gap by disorder, which is in contradiction to 
experiments. Thus correlated defects are required to explain the persistence of the energy gap. 
We have also studied the IPR for various system sizes and find that states in the flat region 
and the upper part of the slope are extended, whereas states in the band tail are localized, as 
expected. Figure ^| shows that the states at the Fermi energy become more and more extended 
with annealing, since the disorder potential decreases. 

In Fig. |] the IPR for fully annealed configurations is plotted as a function of Mn 
concentration x. The number of holes per Mn, p, has been chosen in accordance with 
experiments [[3j]. The change in the IPR-curves with x is rather weak due to the strong ionic 
screening. The main effect is the shift of the Fermi energy due to the varying hole density, 
which drives the Fermi edge into localized states below about 1% Mn. This qualitatively 
explains the metal-insulator transition observed in (Ga,Mn)As [[3j]. 
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Figure 4. IPR as a function of energy for equilibrated configurations as a function of Mn 
concentration x with the number of holes per Mn, p, chosen in accordance with experiments 
[0]. The Fermi energies are also indicated. 



4. Ferromagnetism 

Finally, we turn to the magnetic properties. The bulk magnetisation is composed of 
contributions of the local Mn spins and of the holes. Since the exchange interaction between 
valence-band holes and local Mn spins is antiferromagnetic, these contributions are oriented 
oppositely. However, the hole contribution turns out to be very small. We study the 
magnetisation within a selfconsitent mean-field approximation, which we briefly outline in 
the following. 

The starting point is the Hamiltonian of the coupled valence-band holes and Mn spins, 
H = £ & 4 a c na - J pd £s, ■ S ; , (4) 

no i 

where S ; - are Mn spins (5 = 5/2) and 

S .= £ 4 ct ^(r.)^^,(r.) Vct , (5) 

nOn'o' 

are hole spin polarisations at the Mn sites. £ n and \\f n are the hole eigenenergies and 
eigenfunctions, respectively, for vanishing magnetic interaction, 7 pd = 0. These have been 
obtained in the previous section — note that H contains the full Coulomb disorder potential. 
£ n includes the hole chemical potential. 
We rewrite Eq. @) as 

H = ZZn<toC*o -yE [(^ + S ; ) 2 - (s ; - S,) 2 ] (6) 

no 4 j 

and employ the imaginary-time functional integral. The two spin-squared terms are decoupled 
using Hubbard-Stratonovich transformations with auxilliary fields coupling to s- ± S-. The 
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auxilliary fields are then transformed to new fields h = h+ + ih. , h'- = h+ — z'h. . This leads 

J i iii i i 

to the Lagrangian 

L = Y, + $n)c n o + 4 COS 0, + £ ^ - £(h f • S . + ■ S ; ) . (7) 
na i i ''pd 

To obtain the mean-field theory we employ a saddle-point approximation. We assume h- and 
to be constant in time but not in space, since we wish to retain the disorder effects. The 
holes and Mn spins are then easily integrated out. For the Mn spins this is best done in the 
Hamiltonian formalism. This leads to the grand potential 

with the mean- field hole Green function 

i 

The saddle-point equations are obtained by setting the derivatives of f5Q. with respect to h ; 
and to zero, 

11, = /^ = 7^^^), (10) 

i 

h ! = V m ^-Tl E ^^(^^(^^^(RJ (ID 

" i(0 nan' a' 

with the usual Brillouin function B s (x). To actually evaluate Eq. ( |TTJ ) we diagonalize the 
hole sector of the Hamiltonian or, equivalently, the inverse Green function G l , writing 
G 1 = &g~ l &*, where = — z'co + | nCT is diagonal. We assume the the mean-field Mn spin 
polarisation M ; . to be collinear and choose the z axis along M-. Then G 1 is symmetric and 
6 is orthogonal. Also, G 1 is then already diagonal in spin space so that we can diagonalize 
separately in each spin sector. We denote the components of 6 for spin sector o by ^ m i. a - 

The inversion of G is now trivial, G = &g&* . The frequency sum can finally be performed, 
since does not depend on frequency. Introducing the new eigenfunctions 

«Wiy = EwW'W> (12) 

n> 

Equation < JTT| ) then obtains the form 

h / = J^LMLo) Y*no(*i) Vnofri). (13) 

nd L 

The mean-field solution is obtained by iterating Eqs. ( jTOD and (JT3|). 

Figure |^ shows the magnetisation curves for Mn and hole spins, again for x = 0.05 
and p = 0.3 for random and equilibrium configurations. The hole-spin contribution to the 
total magnetisation is found to be very small. The shape of the Mn magnetisation curve 
becomes more "normal" and concave from below with annealing |JJ], except for a tail at higher 
temperatures, in good agreement with experimental annealing studies [|2|]. The anomalous 
shape for random defects is due to the localisation tendency of the holes, which leads to a 
shorter-range effective Mn spin interaction. This, together with the random positions of spins, 
leads to a broad distribution of effective fields acting on these spins. For higher temperatures, 
only a few spins with strong interactions carry most of the magnetisation, leading to the drawn 
out tail for random defects. The noise apparent on the curves is due to there being several 
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Figure 5. Magnetisation of Mn spins (up) and hole spins (down) as functions of temperature 
forx = 0.05 and p = 0.3 for random, partially annealed, and equilibrated configurations. The 
dotted curves give the standard deviations AM of Mn spins. For comparison, the long-dashed 
curves show the magnetisations obtained from a theory that totally neglects disorder. 

saddle-points with similar free energy and total magnetisation. We also plot the standard 
deviations AM of the Mn spin polarisations. The standard deviation becomes comparable with 
the average Mn spin polarisation at intermediate temperatures below the Curie temperature 
T c . This shows that the reduction of the total magnetisation with increasing temperature is 
initially dominated by disordering of large local moments. Only at higher temperatures the 
local moments also decrease. 

The Curie temperatures are significantly larger than what one obtains from a virtual 
crystal approximation (VCA), neglecting disorder completely, see Fig. |5|. It is an important 
open question whether this is an artifact of mean-field theory or indeed a physical result. 
However, the tails of the magnetisation curves with small average Mn spin polarisation but 
large standard deviation are dominated by a few anomalously strongly coupled spins. These 
are overemphasized by mean-field theory (remember that mean-field theory gives an incorrect 
finite polarisation even for two coupled spins). Thermal fluctuations should be effective in 
destroying the long-range order in this case, thus reducing T c . 

Finally, Fig. || shows magnetisation curves for equilibrated configurations at various Mn 
concentrations with the number of holes chosen as above. We see that the curves for smaller 
Mn and hole concentrations become more similar to the random limit in Fig. |5], which is easily 
understood since the holes become more localized. 

5. Conclusions 

We have shown that strong Coulomb interactions in highly compensated diluted ferromagnetic 
semiconductors naturally lead to strong correlations in the positions of charged defects. We 
also found that these correlations are essential for the description of transport properties, in 
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Figure 6. Magnetisation of Mn spins (up) and hole spins (down) as functions of temperature 
for various Mn and hole concentrations. 



particular the metal-insulator transition, the persistence of the energy gap, and the shape of the 
magnetisation curves. It is a pleasure to acknowledge discussions with P. J. Jensen, J. Konig, 
U. Krey, W. Nolting, M. E. Raikh, F. Schafer, and J. Schliemann. 
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